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ABSTRACT 


An  approximate  solution  to  the  Navier  equations  of  the  three- 
dimensional  theory  of  elasticity  for  an  ex  1  symmetric  orthotropic  circular 
cylinder  subjected  to  internal  end  external  pressure,  axial  loads,  and 
closely  spaced  periodic  radial  loads  is  introduced.  Numerical  comparison 
with  the  exact  solution  for  a  transversely  isotropic  cylinder  subjected 
to  periodic  band  loads  shows  that  very  good  accuracy  is  obtainable. 

When  the  results  of  this  approximate  solution  are  compared  with 

ii 

previously  obtained  results  of  a  Flugge-type  shell  solution  of  a  ring- 
reinforced  orthotropic  cylinder,  it  is  found  that  the  shell  theory  gives 
a  fairly  accurate  representation  of  the  deformation  except  in  the  neighbor¬ 
hood  of  discontinuous  loads.  The  addition  of  transverse  shear  deformations 
does  not  improve  the  accuracy  of  the  shell  solution. 

When  Hill's  orthotropic  yield  criterion  is  applied,  it  is  found 
that  yielding  could  begin  rather  early  at  the  inner  surface  of  the  shell 
adjacent  to  the  frame.  It  is  noted  that  the  transverse  shearing  stress 
has  no  great  effect  on  the  initial  yield  pressure. 
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LIST  OF  SYMBOLS 

elastic  coefficients  in  stress-strain  relations, 
Eqs.  (3*3b) 

arbitrary  constants  of  integration  in  solution  to 
generalized  plane  strain  problem,  Eqs.  (5.8a,b) 

band  load,  Eq.  (2.1a) 

axially  varying  portion  of  band  load  (load 
for  first  problem),  Eq.  (2.2) 

elastic  constants  in  genejralized  plane  strain 
elastic  .law,  Eqs.  (5>2c) 

constants,  Eqs.  (4.13) 

half  width  of  band  load.  Fig.  1 

nondlmens ional  elastic  constants  appearing  in 
Navier  equations,  Eqs.  (3 -5c) 

radial,  circumferential,  ial  and  transverse 
shear  elastic  moduli,  respectively,  Eqs.  (3>2) 

elastic  modulus  of  an  isotropic  material,  Eq.  (3*7) 

exponential  functions  of  shell  geometry  and 
Fourier  index,  Eqs.  (4.13) 

arbitrary  constants  of  integration  of  asymptotic 
solution  to  first  shell  problem,  Eqs.  (4.12) 

nondimens ional  elastic  constants  relating  asymp¬ 
totic  displacement  solutions,  Eqs.  (4. 12g) 
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nondimensional  Fourier  coefficients  of  externally 
applied  loads  and  radial  band  load,  respectively, 
Eqs.  (4.13)  and  (2.1b) 


constants  appearing  in  plane-strain  solutions, 
Eqs.  (5.8) 

center-to-center  distance  between  successive  band 
loads,  Fig.  1 


s  .  $  .  s  ,  $ 

r  0  '  z’  rz„ 
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respective  radtal,  circumferential,  axiel  and 
transverse  shear  stress  coefficients  in  tri¬ 
gonometric  series 


nondimensional  axial  and  radial  displacements, 
respectively,  Eqs.  (4.4) 


x,  e,  z,  s 


radial,  circumferential,  axial  and  transverse 
shear  yield  stresses,  respectively,  Eq.  (6.1) 


radius  of  datum  surface  (usually  mean  radius), 
Eqs.  (4.2) 
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nondimensional  elastic  coefficients,  Eqs.  (4.8e) 


V  V  ez 


classical  radial,  circumferential  and  axial 
strain  parameters,  respectively,  Eqs.  (3.1) 


e 


Naperian  logarithmic  constant,  Eq.  (4.2c) 


plastic  potential  function  of  the  stresses, 

Eq.  (6.1) 

radially  varying  functions  in  the  solution 

for  the  displacements  in  the  first  problem,  l 

Eqs.  (4.4) 

asymptotic  approximations  to  fn  and  g^j, 
respectively,  Eqs.  (4.10) 

»  mTa/4,  Eq.  (4.4c) 

Fourier  series  index  integer,  Eqs.  (2.1) 

uniform  external  and  internal  pressures, 
respectively,  Fig.  1 

magnitude  of  band  load,  Fig.  1 

inner  lateral  surface  pressure  for  plane 
strain  problem,  Eq.  (5.7b) 

pressure  at  which  yielding  begins  at  the 
position  (r,z),  Eq.  (6.2) 

radial  coordinate,  Fig.  1 

internal  and  external  surface  radii,  respectively, 

Fig.  1 

axial  and  outward  radial  displacements,  respectively. 

Fig.  I 

nond i mens  Iona 1  radial  and  axiai  transformation 
coordinates,  Eqs.  (4.2) 


values  of  x  at  outer  and  inner  shell  surfaces, 
respectively,  Eqs.  (4.2) 

axial  coordinate,  Fig.  1 

roots  of  characteristic  equation  of  asymptotic 
Navier  equations,  Eqs.  (4.12) 

constants  appearing  in  generalized  plane  strain 
solution,  Eqs.  (3-6) 

normalized  exponents  appearing  in  asymptotic 
solution  to  first  problem,  Eqs.  (4.12) 

externa 1-to-internal  radius  ratio,  Eq.  (5-8i) 

classical  shearing  strains,  Eqs.  CM) 

-t/awr,  Eq.  (4.7c) 

cylindrical  coordinate,  Eqs.  (3-1) 

nondimens  Iona!  exponent  appearing  in  ortho¬ 
tropic  generalized  plane  strain  solution,  Eq. 

(3-6f) 

Poisson  ratios,  Eqs.  (3*2a)  and  (3*7) 

nondimensional  radial  transformation  coordinate, 
Eq.  (4.6) 
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values  of  5  at  outer  and  inner  shell  surfaces, 
respectively 

racial,  circumferential,  axial  normal  stresses 
and  transverse  shear  stresses,  respectively 
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1. 

1.  INTRODUCTION 

To  date,  Investigators  of  pressurized  ring-reinforced  circular 
cylindrical  shells  have  employed  shell  theory  in  order  to  describe  the 
behavior  (see,  e.g.,  Ref.  I).  For  thin-walled,  isotropic  shells,  in  which 
the  collapse  mechanism  is  primarily  a  buckling  instability,  the  results  of 
shell  theory  are  generally  quite  satisfactory.  However,  in  short,  thick- 
walled,  filament  wound  composite  cylinders,  the  collapse  mechanism  is  quite 
complicated  (as  yet  no  single  description  has  been  agreed  upon).  Examina¬ 
tion  of  the  ruptured  composite  test  models  indicate  that  failure  is  prob¬ 
ably  due  to  a  transverse  shear  build  up  at  the  frame  (Ref.  2),  a  phenomena 
not  predicted  by  shell  theory. 

Because  of  the  complicated  failure  mechanism  of  composite  shells, 
the  need  is  obvious  for  an  analysis  which  would  yield  a  more  accurate  des¬ 
cription  of  the  stress  distrioution  through  the  thickness  than  is  presently 
available  through  existing  shell  theories.  The  work  presented  by  Klosner 
and  Levine  (Ref.  3)  on  isotropic  shells,  and  later  extended  by  them  (Ref.  4) 
to  include  transversal ly  Isotropic  shells,  consisted  of  an  elasticity  analy¬ 
sis  in  which  the  stress  functions  suggested  by  Leknitskii  (Ref.  5)  were  used 
to  satisfy  exactly  the  classical  equations  of  elasticity.  Their  results 
demonstrated  clearly  that,  although  the  commonly  used  Donnell  shell  theory 
gave  excellent  results  throughout  most  of  the  shell,  it  could  not  accurately 
predict  the  stress  distributions  at  the  frame.  However,  the  Leknitskii 
stress  functions  are  useful  only  for  the  case  of  a  transversal ly  isotropic 
(or  isotropic)  material.  The  cyl indr  Seal ly  orthotropic  nature  of  the  filament 
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wound  composite  shells  can  only  be  accounted  for  by  seeking  a  solution  to  the 
more  complicated  equations  of  elasticity  for  such  a  material.  In  the  present 
study  an  approximate  solution  for  the  displacements  is  used  to  solve  the 
axisymmetric  Navier  equations  of  equilibrium  for  an  infinitely  long,  ortho¬ 
tropic,  hollow  cylinder  under  external  uniform  pressure  and  closely  spaced 
periodically  varying  internal  band  loads.  It  is  indicated  how  more  accuracy 
may  be  obtained  by  use  of  either  a  perturbation  or  iteration  method. 

it  is  found  that  the  asympto  ’c  solution  is  quite  accurate  and 
that  the  shell  solution  presented  in  Ref.  I  is  in  very  good  agreement  with 
the  asymptotic  solution,  except  in  the  vicinity  of  the  frame.  When  trans¬ 
verse  shear  deformations  are  considered  in  the  shell  theory  no  better  agree¬ 
ment  is  obtained. 


^  , 
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2.  FORMULATION 


The  shell  of  infinite  length  is  subjected  to  an  external  hydrostatic 
pressure  load  and  internal,  axisymmetr ic,  periodically  spaced  band  loads. 
These  band  loads  are  described  as  (e.g.,  see  Ref .  4) 


GO 

B  =  p  [  £  H  cos  -  A] 
C  C  n=l  n  1 


(2-la) 


wi  th 


H„  “  m  H,"+'sin 


and 


(2.1b) 


A  =  D/l 


(2-lc) 


where  pc  is  the  magnitude  of  the  band  load,  D  is  half  the  width  of  the  band 
load,  t  is  half  the  distance  between  the  centers  of  two  successive  band 
louJ.ij  and  z  is  the  axial  coordinate  (see  Fig.  1). 

When  the  classical  linear  equations  of  three-dimensional  elasticity 
are  used  to  describe  the  shell  deformation,  the  solution  may  be  taken  as  the 
superposition  of  the  following  two  solutions: 

1.  The  solution  of  the  shell  subjected  to  an  internal  load  des¬ 
cribed  by  [see  Eqs.  (2.1)] 


B  .  =  p  EH  cos 
cl  rc  .  n 
n=l 


DEL 

l 


(2.2) 
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corresponding  to  which  the  displacements  can  be  taken  as  trigonometric  series  I 

in  z.  In  the  ensuing  analysis  the  classical  elasticity  equations  are  used  ' 

to  find  the  coefficients  in  the  latter  series  as  functions  of  the  thickness 
coordinate  r. 

2.  A  generalized  plane  strain  solution  which  includes  uniform 

lateral  pressure  loadings  on  the  inner  and  outer  surfaces  [the  inner  surface 

load  must,  of  course,  include  the  term  ■  p  i  appearing  in  Eq.  (2.1a)]  as 

c 

well  as  a  constant  axial  force. 


3-  BASIC  EQUATIONS 


The  classical  strain-displacement  equations  used  to  describe  the 
axisymmetric  deformations  of  a  circular  cylindrical  shell  are 


e  =  w 


e  =  u 

2  ,2 


7r2SV  +  w,2 


eQ  =  w/r 


yr0=  0 


^2=° 


(3- la, b) 


(3-lc,d) 


(3- le,f) 


in  which  a  comma  followed  by  a  subscripted  variable  denotes  differentiation 
with  respect  to  that  variable,  r,  Q  and  2  are  the  radial,  circumferential 
and  axial  coordinates,  respectively,*  er,  e0  and  ez  are  the  corresponding 
normal  strains;  yfQ}  y ^  and  yQz  are  the  shear  strains,  w  and  u  are  the 
displacements  in  the  radial  and  axial  directions,  respectively.  It  should 
be  noted  that  both  r  and  w  are  taken  positive  outward  (Fig.  1). 

The  generall2ed  Hooke's  law  for  a  cyl indrical ly  orthotropic, 
homogeneous  material  is  taken  as 


’  vre 


■V  '  1 


*Vrz\  ar 


"V02  4ae 


(3-2a) 
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T  =  G7 

rz  'rz 


n  re 


=  0 


Tez  =  0 


(3.2b,  c,d) 


where 


) 


} 


(3.2e) 


follow  from  symmetry  of  the  stress  and  strain  tensors.  {  }  denotes  a  column 

matrix  and  (  )  denotes  a  square  matrix,  a  ,  oa  and  0  are  the  normal  stresses 

r  y  z 

in  the  radial,  circumferential  and  axial  directions,  respectively;  rrz  is  the 
transverse  shear  stress.  E^,  Eg  and  Ez  are  the  elastic  moduli  in  the  radial, 
circumferential  and  axial  directions;  G  is  the  transverse  shear  modulus 
and  vr0,  v0r,  vr2,  vzr,  vflz  and  vz0  are  Poisson  ratios.  It  can  be  seen 
that  Eqs.  (3-2)  contain  7  independent  elastic  constants. 

The  inverse  relation  is  given  by 


where 


(3.3a) 


-  0-v.aVo,)(e./n) 


20  ez' 


A22  *  <'-v«vzr)f6«/N> 


A33  =  (|-vr9v9r)<£i/H> 


(continued  on  next  page) 
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*12  -  A2I  *  *  (vrfl  +  vrzvz0)(E8/H)  <  *13  “  *31  ’  '  <Vrz  *  Wez*  !EzA!) 


*23  *■  *32  •  ‘  <v8z  +  v8rvr2>  (Ez/N) 


"  *  1  -  vz8v0z  -  v8r(vr8  +  vrzvz9)  '  vzr(vr8v8z  +  vrz> 


(33b) 


The  axisymmetric  equilibrium  equations  are 


(ra  )  -  au  +  rr  .  =0 

'  r ,  r  6  rz,z 


(Ma) 


r<r  +  (tt  )  =0 

z,z  '  zr  ,r 


(3.4b) 


Equations  (3.i),  (3-3)  and  (J.4)  are  combined  to  yield  the  well- 
known  Navier  equations  for  an  orthotropic  cylinder  undergoing  axisymmetric 
deformations  (Ref.  5)- 


r\rr  +  rw,r  *  “/".zz  *  °2“  +  °3rSrz  +  Vu,z  '  0 


(3-5a) 


and 


2  2  2 
D«r  w  -  Dy-rw  +  D7r  u  +  r  u  +  ru  =  0 
y  )  o  ,z  /  ; zz  y  x r  >  r 


(3.5b) 


where 
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0,  =  G/A,  j  ,  Dk  =  (A,„  -  AlX)/Aj ,  ,  D,  =  A„/G 


l"23  n]3,,n\\ 


7  "  "33' 


°2  =  A22/A 1 1  '  V  (G  *  V/G 


°3  =  (G  *  A13>/A1 1  »  d6=(a32-g)/g 


(3.5c) 


It  is  noted  that  only  five  of  the  constants  Dj,  *•*,  are  independent,  since 


0j  =  D,Dj  and  04  -(Dj  +  Og)o, 


(3.5d) 


Equations  (3>5a)  and  (3.5b)  are  a  set  of  second  order,  linear,  homogeneous 
partial  differential  equations  with  variable  coefficients. 

If  the  material  is  transversely  isotropic, 


Er  -  Ee  ’  vrS  ■  V  ’  v”  *  v«.  ’  v~  *  v- 


Ai*  =  A*,  =  A»o  ,  A, ,  =  A, 


rz  r0z 


zr  rz  9 


‘13  "  31  "  "32 


'll  ”  22 


02  =  1  ,  D5  =  -  °6  >  °4  =  0  (5.6) 


It  follows  that  the  number  of  independent  elastic  constants  reduces  to 
five. 

If  the  material  is  isotropic, 


Er  ■  Ee  ■  Ez  ■  E  '  vr9  V  *  vrz  “  V  ”  vzr  '  *  v 


(continued  on  next  page) 
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4.  AXIALLY  VARYING  LOADS 

The  first  problem  to  be  investigated  is  that  of  an  infinitely  long 
circular  cylinder  with  an  internal  load  described  by  Eq.  (2.2).  For  such  a 
loading  the  solution  is  any  solution  to  the  Navier  equations,  Eqs.  (3-5); 
which  satisfy  the  conditions  that: 

on  the  outer  surface  of  the  shell 


a 

r 


0 


T  rz  *  0 


and  on  the  inner  shell  surface 


mzz. 


a  =  p  E  H  cos  . 
r  Kc  ,  n  L 
n=  I 


Trz  =  0 


(4.1a) 

(4.1b) 


(4.1c) 

(4. Id) 


The  major  difficulties  in  finding  such  a  solution  are  caused  by  the  fact 
that  the  coefficients  in  the  Navier  equations  are  variable  in  r.  In  order 
to  reduce  the  number  of  variable  coefficients  encountered  (and  at  the  same 
time  nondimensional ize  the  equations)  the  following  well-knovai  transformation 
is  introduced: 


x  =  (r/a) 


(4.2a) 


' 

j 

i  ■ 


\ 


f 

*  i 


n. 


z/a 


(4.2b) 


with  the  Inverse  relation 


r  =  ae 


(4.2c) 


z  =  ay 


(4. 2d) 


in  which  a  Is  the  radius  to  any  selected  datum  surface  such  that  r.  <  a  <  r  . 

i  —  “  o 

Hence,  for  relatively  thin-walled  cylinders  x  =  (r/a)  -  1,  which  is  small 
compared  to  unity.  The  Navier  equations  now  become 


W,*x  +  D!e2Xw,yy  -  °2U  +  °/U,xy  +  \<>\y  *  0 


(4.3a) 


D-eXw  -  D/-eXw  +  D^e^u  +  u  =0 
5  ,xy  6  ,y  7  ,yy  ,xx 


(4.3b) 


The  separation  of  variables  method  is  used  to  find  a  solution. 

The  assumed  from  of  the  displacements  must  reflect  the  symmetry  of  the  load¬ 
ing  (Fig.  1).  Thus, 


A.  .W  oo 

c  n=  l 


(4.4a) 


u  =  5r~=  \ 


(4.4b) 


where 


k  =2P 

n  1 


(4.4c) 


and  W  and  U  are  nondimensional  displacements.  Hence,  the  Navier  equations 
yield 


g  -  (kZD.eZx  +  Djg  +  k  D,eXf  +  k  D.eXf  =  0 
sn,xx  v  1  2  3n  n  3  n,x  n  4  n 


2x, 


m  k  D  c  o  k  Dse  q  •  k  0  c  f  +  f  ~  q 
n  5  n,x  n  6  3n  n  7  n  n,xx 


(4-5a) 

(4. 5b) 


A  significant  insight  is  gained  into  the  nature  of  the  functions  gn  and  f^ 
upon  the  stretching  of  the  x-coordinate  through  the  introduction  of  the  trans 
formation 


|  =  kRx  =  (mra/^)x 


(4.6) 


The  Navier  equations  finally  yield 
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g  £  -  (D.e  +6  Djg  +  D,e  "  f  .  +  6  Dke  f  =  0 

sn,||  '  1  n  2,3n  3  n,|  n  4  n 


(4.7a) 


Equations  (4.7a, b)  are  second  order,  linear,  simultaneous, 
ordinary  differential  equations  with  variable  coefficients-  These 
equations  can  be  uncoupled  with  the  aid  of  integrating  factors-  How¬ 
ever,  the  resulting  relations  have  complicated  variable  coefficients 
and  are  singular  for  the  not  uncommon  case  o>'  a  transversely  isotropic 
cylinder.  The  solution  obtained  for  transversely  isotropic  cylinders 
in  Ref.  4  (which,  as  mentioned  previously,  is  formulated  in  terms  of 
stresses  and  solved  with  the  Lehknitskii  functions)  is  in  terms  of 
Bessel  functions-  Accordingly,  if  Eqs.  (3*6)  are  introduced  into 
Eqs.  (4.7a, b)  and  the  resulting  relations  are  uncoupled,  the  displacement 
solution  is  also  in  terms  of  Bessel  functions,  which  correspond  to  those 
found  in  Ref.  4. 

The  success  of  Klosner  and  Levine  (Ref. 4)  in  using  Bessel 
functions  for  the  transversely  isotropic  case  tempts  one  to  use  a  Frobenius 
or  other  type  of  power  series  solution  for  the  more  complicated  orthotropic 
equations.  However,  §  grows  in  proportion  to  n  [Eq.  (4.6)]  and  unless  (as 
In  classical  shell  theory)  the  thickness-to-length  ratio  is  very  small 
and  the  load  converges  quickly  (i.e.,  is  continuous  and  varies  slowly  with 
y)  <•  becomes  large  as  n  increases.  The  large  values  of  kn  (and  hence  of  |) 
which  evolve  from  the  short  shells  of  interest  here  make  a  power  series 
solution  impractical  in  that  many  terms  are  necessary  for  convergence. 

Also,  such  solutions  may  give  rise  to  small -difference  terms,  unless 
appropriate  asymptotic  expansions  can  be  obtained. 
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On  the  other  hand,  when  kn  is  large,  is  small  [Eq.  (4.7c)]. 
This  suggests  that  a  perturbation  or  iteration  scheme  may  yield  a  reason¬ 
able  approximate  solution.  Any  solution  which  capitalizes  on  the  size  of 
6^  might  also  be  thought  of  as  an  asymptotic  expansion. 

Once  and  f  have  been  found,  the  stresses  can  be  determined  by 
^n  n 

combining  Eqs.  (4.4)  and  (3-3) •  Hence, 


0  od  ,  -6  £ 

1  r  *  rr 


n*  r  [e  "  <9n,s  -  al26nSn)  -  a|jf„)cOSkny 


-  -  Z  l-e’6nE''al29n,S  ‘  a22S9„>  ‘  a23fn)c0skny 
pc  n=l  n 


°z  00  1 

— -  =  E  ?— (-  e  (a,,g  ,  +  aoX5  g  )  +  ax*f  ]cosk  y 

Pc  n_i  5n  I33n, £  23  n3n  33  nJ  n7 


Trz  ”  1  r  8n^f  ,  -  g  jsink  y 
- —  =  E  r”  le  n, |  an  n7 


P  ”,  6 
*c  n=l  n 


whe  re 


(4.8a) 


(4.8b) 


(4.8c) 


(4.8d) 


3 12  "  A 1 2/A 1 1  7  al3  "  °1  ■  °3  7  a22  “  °2 

a23  =  °4  +  Di  -  d3  7  a33  =  D7°l  ^‘8e) 
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The  boundary  conditions  are  satisfied  by  [see  Eqs.  (4.1)  and 
Eqs.  (4.8a,d)] 


a  °(9n,i  ‘  *l28n3n)  "  a15f„  =  0  (l*-9a) 


(4.9b) 


at  the  shell  outer  surface  x  =  xq,  |  =  &  >  0 


and 


-x, 


e  f  -  g  =  0 

M  n 


-x , 


e  (gn,5  ’  alA9n>  *  aljfn  =  Hn 


at  the  shell  inner  surface  x  =  x^  |  =  |j  <  0 


(4.9c) 


(4  9d) 


If  a  perturbation  solution  is  used  in  which  the  functions  gn  and 

f  are  considered  to  be  perturbed  about  about  the  solution  of  Eqs.  (4.7) 

with  6  set  to  zero,  g  and  f  may  be  expressed  as 
n  n  n  r 


CO 


9" -  ,1 

9  •  (5 )& 
an  i  n 

(4. 10a) 

CO 

f  =  I 
n  i=l 

f  •  (1)6* 

ni  '  n 

(4. 10b) 
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For  i  =  1  (zerot^1  perturbation)  the  differential  equations  are 


3nl,§?  "  D 1 9n  1  +  D3fnM  ° 


(4.11a) 


a"  d 


Dcg  ,  ,  +  f  ,  “  D-f  ,  =  0 

53n1,£  nl,5§  7  nl 


(4.11b) 


Equations  (4.11)^  in  which  6n  =  -k/nrra  =  0,  correspond  to  the  equa¬ 
tions  of  a  doubly  infinite  flat  plate  (for  which  a  -♦  ®)  with  a  thickness 
coordinate  £  =  ax.  The  solution  is  written  as 


4  0. 

q  ,  =  £  G^.e  1 

i  =  l 

4  p" 

f  ,=  E  H..G"  e  1 
n  1  1 1  1 1 


(4.12a) 


(4.12b) 


where 


P",2  *  al,2(5  -  S0) 


(4.12c) 


@3,4=  -a2,l(5  '  ?l> 


(4. I 2d) 


and  a,  and  are  those  roots  of  the  characteristic  equation 


a)  +  (D305  -  0?  -  D,)af  +  d,07  =  0 


(4. 12e) 


which  are  defined  by 


a 


l,2=-l(D3D5“D7-01)±i[(D3D5-07-DJ)  -  40,07] 


1/2 


(4.12f) 
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Furthermore, 


a,  -D,  D_at 


Hu  -  -  (-ten  =  <  2  > 

1  3  i  ai  “D7 


a3,4  =  ‘a2,l 


(4.12g) 


(4. 12h) 


For  the  specific  material  properties  considered  subsequently, 

«l  is  real;  moreover  Oj  2  >  0  and  a.+j  %  a. .  In  general,  the  nature  of  a. 
is  dependent  only  upon  the  elastic  properties  of  the  shell  material.  For 
example,  if  the  shell  were  isotropic,  equal  roots  for  the  characteristic 
equation  would  be  found  (i.e.,  Gj  ^  =  “  ')•  It  should  also  be  noteo  that 
the  roots  of  the  characteristic  equation  do  not  change  with  n.  The  exponents 
p"  were  selected  to  modify  G^.  in  order  to  avoid  numbers  of  excessive  magnitude 
in  the  computations.  The  arbitrary  constants  G^.  are  obtained  by  satisfying 
the  boundary  condition  Eqs.  (4.9)  which,  when  combined  with  Eqs.  (4.12) 
yields  the  matrix  equation 


(cfjUGf,}  -  {H? }  ,  i,j  =  1, 2,3,4 


(4.13a) 


where 


CU  "I*  °<VSnal2>  -al3Hn)Ell 


c?.  -  (e  °a.H,.  -  l)E" 
2i  '  ill  'll 


(4. 13b) 


(4. 13c) 
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C3i  ”  ^  aiHli-1^E2i 


(4.13d) 


C4i  =  1(ai"a126n^  "  a 1 3H 1 i J E2 i 


(4. 13e) 


_n  _  pn  _  _n  _  Fn  . 

1 1  12  fc23  24  ~  1 


(4.130 


Pn  Fn  WU  . 

EI3  “  E22  =  e 


cn  _  cn 

El4  E21  “  e 


“iMo*  aikn^xrxo^ 


(4. I3g) 


(4. 13h) 


For  the  loading  of  interest  here  (see  Eqs.  (4.9)  and  (4.13)] 


Hi  =  Hj 


h3  =  0 


(**•  130 


H  =  H 


(4. 13j) 


Of  course,  these  conditions  may  be  easily  altered  to  include  other 
loading  conditions}  a  more  detailed  discussion  on  this  and  a  means  to  inves¬ 
tigate  layered  shells  is  offered  in  the  Appendix. 

Once  the  constants  .  have  been  determined  from  Eq.  (4.13a),  the 
stresses  may  be  obtained  by  substituting  Eqs.  (4.12)  and  (4.10)  into  Eqs. 
(4.8)}  this  yields 


CF  or>  4  (3  .  oo 

~  =  £  E  (e  X^a. -aj^J-a^d]  j  je  'cosk^v  =  £  cosk^y 


c  n=l  i  =  1 


n'  ,  r  n‘ 
n=  1  n 


(4.14a) 


19. 


£■-  r,  A  1-e’*(*U0,|-*22Wli1Gi/l“skny=  r,  Secoskny 


c  n=l  1=1 


n=l  n 


„*  .f|['e'X(a13ai+a23Sn)  +  a33HU)GlieS'COSkny  *  i,S2„C°5kny 


n=l  n 


T  OB  If  £ .  » 

—  =  D,  E  E  (e"Xa.H,  j-ljG^.e  'sink  y  =  E  S  sinky 
Kc  n=l  i= 1  n=1  n 


(4. l4b) 


(4.14c) 


(4. I4d) 


If  more  accuracy  is  desired,  the  first  perturbation  can  be  applied 
[  i  -  2  In  Eqs.  (4.10)  and  (4-7)].  This  yields 


9n2,U  "  °l9n2  +  °3fn2,|  "  *2Dl9n l"D3fn l,!*1  “  Vnl 
*D59n2,§  +  fn2,K  "  Vn2  “  (°591,§  +  2°7flH  "  °69ni 


(4. 15a) 


(4.15b) 


An  iteration  technique  may  also  be  used.  The  first  iteration  con¬ 
sists  of  Eqs.  (4. 12),  which  corresponds  to  &n  =  0.  Any  further  iteration 
(gni,  fn.)  may  be  obtained  by  successively  solving 

26  £  -  6  £  6  £ 

9  •  t * “D 1 9  -+D_f  .  =  (D.  (e  n  -l)+6  D,]g  .  ,+D_(l-e  n  )f  .  ,  -e  n  6  D. f  .  , 

3m,££  1  ni  3ni,|  P  n  2  am-l  3  '  m-l,£  n  4  m-l 

(4. 16a) 


6  §  5J  26  £ 

-D_g  .  .+f  .  e, -D„f  .=DR(e  n  -l)g  .  .  .-5  D^e  n  g  .  ,+D„(e  -l)f  .  . 
5  ni,£  m,££  7  ni  5  m-l,|  n  6  3m-l  7  n  1  - 1 


(4. 16b) 


20. 


where  the  right  hand  sides  are  obtained  by  inserting  results  of  prior  iterations. 

in  either  the  perturbation  or  iteration  technique  the  stress  components 
corresponding  to  perturbations  or  iterations  beyond  the  zeroth  must  vanish 
on  the  boundaries.  In  either  approach,  a  correction  term  is  only  necessary 
in  the  early  terms  of  the  trigonometric  series  given  in  Eqs.  (4.4),  (i.e., 
when  n  is  relatively  small),  since  &n  becomes  smaller  as  n  increases  and, 
therefore,  the  approximate  solution  represented  by  Eqs.  (4.12)  becomes  more 
accurate. 

The  selection  of  p1?  to  modify  the  constants  is  indeed  appropriate, 
for  as  n  grows  large  the  coefficients  may  be  closely  approximated  by  simple 
formulae  which  reveal  the  dependent  of  the  trigonometric  series  convergence 
on  location  through  the  shell  thickness.  Hence,  as  n  becomes  large,  kn  be¬ 
comes  large  and  [see  Eqs.  (4.13g)  and  (4.13h)) 

E]3,  E22,  E]V  E2)  -  0  (4.17a) 


Therefore, 


»n  aH  «n  pn  pH  pn  a 

H3'  ul4'  '‘ly  u24>  °31'  1'32’  U4P  u42  "*  u 


and  Eqs.  (4.13)  may  be  replaced  by 


o\ 


0 


0  0 
„n  n 

C33  C34 

Cn  cn 
c43  C44 


Gn 
G11 

rn 

<;av 

G13 


(4. 17b) 

i 


(4.18a) 


V 


Thus,  for  large  n  (n  >  N  say) 


°>>  *  ="a  •  • 


(4.18b) 


_ 

«n  34  n 

u,-.  =  -  ji  jt  j, 


13  ^33^3^: 


(4. 18c) 


'33c44“c43c44 


(4. 18d) 


n  ■«*;(£ -S,)  n  -<*i  (S-S ,)  H 

gn  1  =  *■  "C34e  +  c33e  1  (cn-crt  -cn  cn  ^ 

J  3J  c33  44  C43C34 


n  n  -a,  (S-l.)  H 

fnl  -  [W  -  C33H 1 1 e  1 


'33  1 1 


Sr  =  *'"C34(“a2e  X+ai3H]2^e 


(4.19a) 


(4. 19b) 


n  -x  -a, ($-6.)  H 

+  c33^,e  +.,3H|t)e  K-^q,^) 


.  n  ,  -x  *^2  i ) 

Sen  =  ['C34(a12a2e  +a23H12^e 


(4.19c) 


.  n  -x  ""^2(£"£i) 

Szn  ~  “C3lf ] 3^2®  "a33H]2^e 


(4. I9d) 


n  ,  -x  -a. (£-6.)  H 

+  c33(al3c,l'  "a33HH)e 


(4. I9e) 


22. 


Srz  =  t-c"4(a2H12e""-l)e 
n  J 


-x 


-a2^^i)  n  ,  -x  ,v  ““l^'M 
+C23^aiHi]8  -,'e 


n 

1(^w 


(^•190 


For  the  problems  under  investigation  in  this  report  the  stresses 
converge  as 

-a-U-l,) 

e  H  cosk  y 

n  n 1 

The  convergence  is  fastest  on  the  shell  outer  surface  (£  =  6  >  0)  and 

slowest  on  the  shell  inner  surface  (£  =  |j  <  0),  where  i.  matches  the 
convergence  of  the  applied  load  [Eq.  (2.2)]. 
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5-  SECOND  (GENERALIZED  PLANE  STRAIN)  PROBLEM 


The  second  problem  to  be  considered  is  that  of  an  infinite,  un¬ 
supported  shell  subjected  to  internal  and  external  pressure  and  an  applied 
external  axial  force.  If  the  radial  displacement  is  assumed  to  be  only  a 
function  of  r,  and  if  the  axial  strain  is  assumed  to  be  constant,  the  strain 
displacement  equations  are  given  by  Eqs.  (3-1)  in  which  now 


e  =  u  =  const 
z  >z 


V  =  Vze  =  ver 


Hooke's  law  [Eq.  (3* 2a)]  reduces  to 


er  -  8n°r  •  B12ae  -  Vz 


ee  ■  "  B21°r  +  822°e  '  B23ez 


where 


Bll  "  E  ^'"Vrzvzr^  >  ?!2  B21  E  ^Vr0  +  VrzVz0^ 


B13  Vzr  ’  B22  "  Efl  (!_v0zvz0^  >  B 


and  the  equilibrium  equations  reduce  to 


23  “  Vz0 


P. 

<  t 


(5-la) 


(5-lb) 


(5-2a) 


(5-2b) 


(5-2c) 


24. 


(ror>, r  *  “a  *  0 


(5-3) 


The  governing  differential  equation  is  obtained  from  Eqs.  (3*1) 
and  (5*0  to  (5-3) >  thus, 


rlr(ro  )  J  -  (ro  )  =  (^4^)  re 
r  3  r  ,  r  0^2  r  2 


(5-4) 


Also,  the  displacement  equation  may  be  obtained  by  using  the  plane  strain 
assumptions  in  Eqs.  (3*5);  hence, 


r  w  +  rw  -  D0w  =  -  Di.e  r 
,  rr  ,  r  2  4  z 


(5-5) 


It  should  be  noted  that  no  matter  which  of  Eqs.  (5*4)  or  (5-5) 
is  utilized,  the  previously  mentioned  singularity,  which  occurs  for  either 
a  transversely  isotropic  or  an  isotropic  material,  persists  in  the  generalized 
plane  strain  equations  (i.e.,  when  =  0  and  D2  =  1  or  B2g  =  B^  and 
B  ]  j  =  822).  This  singularity,  which  determines  whether  or  not  the  differen¬ 
tial  equation  is  homogeneous,  is  reflected  in  the  solution.  The  solution 
of  the  present  problem  is 


.  4-1  .  .  -4-1  * 

a  =  A.r  +A„r  -ae 
r  1  2  0  z 


aQ  =  X  (A  l  r*’" '  “a2  r~k~ ' )  ~(x‘oez 


(5.6u) 


(5.6b) 


°z  =  (B13+x'B23)AlrA’"1+(813'XB23)A2r"X"1+[Ez^o(8l3+B23)]ez  (5<6c) 


i  ; 


D,  * 

w  =  (X822-821)A1rX-(XB22+B2j)A2r‘X-(Tf5-)rez 


(5-6d) 


u  =  ze 


(5-6e) 


in  which 


X  =  i  8]]/B22  “  2 


(5.6f) 


(B|3-B23)/(1-Xz)B22  ,  when  X  *  1 


0,  when  X  =  1  (transversely  isotropic  or  isotropic) 


(5-6g) 


%  * 
(W  ‘ 


0j^/(1-D2),  when  X  4  1 


0,  when  X  =  1 


(5-6h) 


The  three  conditions  to  determine  Aj,  A2  and  ez  are  the  boundary 


conditions 


a  =  -  p  at  r  =  r  (  outer  surface), 
r  o  o 


(5.7a) 


Tr  =  -  pj  =  -  (Pj  4.  ^pc)  at  r  =  r]  (inner  surface)  (5-7b) 


and  the  prescribed  axial  force  condition 


f  a  rdr  =  =-  (r.^p.  -  r  Zp  ) 
J  z  2  '  1  *i  o  vo' 


(5-7c) 


I 


26. 


Substitution  of  Eqs.  (5-6)  into  (5-7)  results  in  three  linear  algebraic 
equations  which  may  be  solved  to  yield 


A,  =  Jj/J  ,  A2  =  J2/J  ,  Ez  =  J3/J 


(5-8a,  b,c) 


where 


J|  -  t-8?  ,  t  “*@2  ( y 1  -x- 1 )  - 1 )  r"" 1  J+p0  toysj  ( r 1 (r^  - 1 )  3  j  r . 

f|  (5-8d> 


2  i  \  _  ~X.~  1 .  .  r  „,v *«  "  / .  1  “X 


1-X 


J,=t\  (i-7*-',)<+p)la>|(rul.|)-v3/'l(,2.i))-p  lp|a‘(7wl-i)^3(r2-i)nr.w 

r, 

‘  (5-8e) 


M  . 

J3=  ’ 


I  X.™  I  \  r  ]  /  X’f  1  ,  V  Q  '*  X"  1  /  1  “X»  |  \  1  frt  /  X*f  1  |  \  *'*  /  ]  "X  1  \  1 

-7  )+P,[P,7  (7  -0-P27  ^  -l)]-P0[Pj(7  -0-P2(7  -1)1 

(5-8f) 


J  =  O^|(7W'-l)(r‘X‘-l)^(7l'*-'l)(7X'10^3(72-l)(7‘!l'l-7>-'1)  (5.8g) 


M  1  ,  2  2  , 

M  =  t-  (r.  p.  -  r  p  ) 
p  2  '  1  ri  o  ro' 


(5-8h) 


(5-8i) 


b  -  1l ilifa 

P1  1  +  X 


(5-8j) 
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f(813“X'B23>/(1  ”X)  '  when  x  *  1 


0  ,  when  A,  =  I 


(5-8k) 


a“(B.,  + 
ov  13 


62j)] 


(5-8t) 


The  stresses  and  displacements  of  a  pressurized  cylinder  subjected 
to  periodically  spaced  band  loads  may  now  be  obtained  by  superposing  the 
individually  developed  solutions  [see  Eqs.  (4.4),  (4.8),  and  (5-6)];  this  re¬ 
sults  in 


w 


(XB22-S2  j  )  A ,  rX-  (XB22+52 ) )  (j±-)  ‘“'re.+ 


00  4  B? 

E  6  E  G^.e  'cosk  y 
.  n  .  .  1  i  n 7 


ap  ®  4  0? 

u  =  e  ay  +  I  6  I  H-.G^.e  'sink  y 

z  '  A..  ,  n.  .  1 i  I  i  n' 


(5- 9a) 
(5- 9b) 


0n 
.n  i 


0r=A/'1+A2r'I*,^oVPc  S,  .E.  te_X  (ai  ”a l2^n)~a ljH,  j  ]G|  j e  'cos^y  (5-9c) 

n=l  i=! 


°0  =  X(AlrX"'-A2r"X"1)-<3:oez+Pc  E,  .S,  ['e"X(a12ai‘a225n)"a25Hli]Glie  'C°sknV 

n= I  i=I 


(59d> 


■"z  •  <B.3+XB23)AlrX"'  +  (8l3-^23)A2r'"''  *  ^z 

a  4  0? 

+  Pc  n=l  i=|  i"e"X(al3ai+a235n)  +  d33HlijGlie  'coskny 

®  4  p? 

Trz  a  pc°i  E,  .r.(e"XaiHli'l)Glie  'sinkny 
n=  I  i  =  l 


(5-9e) 


(5-9f) 
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6.  YIELD  CRITERION 


Any  elasticity  solution  is  valid  oniy  up  to  the  pressure  at 
which  yielding  begins  at  some  point  on  the  shell.  It  is  therefore  of 
great  interest  to  attempt  to  find  this  initial  yei Id  pressure. 

R.  Hill,  Ref.  9  >  developed  a  plastic  potential  function  f 
for  the  determination  of  that  combination  of  stresses  for  which  an  ortho¬ 
tropic  material  would  begin  to  yield.  This  function  of  the  stresses  is 
defined  to  be 


f  ■  k  [  <-T  +  -2  -  T>  <V0z,2+(-2  +  T  -  T>  ‘  V°r>2 
2  ©2  z2  x2  9  2  z2  x2  ©2  z  r 


+  (JL+  JL 
T  '  2  2 
X  0 


I  w  v2  tQz  Tzr  Tr0, 

TH"r-"8)  +-r  *~T 

L  R  5  I 


(6.1) 


in  which  X,  0,  Z  are  the  radial,  circumferential  and  axial  yield  stresses, 
respectively,  and  R,  5,  T  are  the  appropriate  shear  yield  stresses-  The 
criterion  for  yielding  to  begin  is  that  f  =  1. 

The  stresses  in  Eqs.  (5-9)  could  be  substituted  to  leave  an  expres¬ 
sion  for  the  pressure  po^(r,z)  for  which  yielding  would  begin  at  the  position 
(r,  z) 


(  \  PQV  /  1  .  1  K  f°e  %2  .  1  1  1  u°z 

Poy(r>2>  ■  p  7)(r'r’  (^+?- j1(r 

1  rc  0  Z  X  rc  rc  Z  X  0  rc 


1  .w>  10x2  ,  'zr  ,-1/2 


+  ^  +  VTH”Tr+"t25 

X  0  Z  pc  pc  P  S 

rc 


(6.2) 
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6.  YIELD  CRITERION 


Any  elasticity  solution  is  valid  only  up  to  the  pressure  at 
which  yie’ding  begins  at  some  point  on  the  shell.  It  is  therefore  of 
great  interest  to  attempt  to  find  this  initial  yeild  pressure. 

R.  Hill,  Ref.  9  ,  developed  a  plastic  potential  function  f 
for  the  determination  of  that  combination  of  stresses  for  which  an  ortho- 
tropic  material  would  begin  to  yield.  This  function  of  the  stresses  is 
defined  to  be 


f  =  2  1  <■?  +  "2  '  T)  <V°2>2+<T  +  T  -  T>  <V°r>2 

2  ez  zz  xz  9  2  l  xz  ez  2  r 


,  /_!_  J.  _I  \  !rl  _  \2  lez  Tzr  Tre, 

+  (  o  +  2  2'  ^°r”a0  +  2  +  2  +  2 

X  ©  Z  r  y  R  S  T 


(6.1) 


in  which  X,  ©,  Z  are  the  radial,  circumferential  and  axial  yield  stresses, 
respectively,  and  R,  S,  T  are  the  appropriate  shear  yield  stresses.  The 
criterion  for  yielding  to  begin  is  that  f  =  1. 

The  stresses  in  Eqs.  (5-9)  couid  be  substituted  to  leave  an  expres¬ 
sion  for  the  pressure  po^(r,z)  for  which  yielding  would  begin  at  the  position 
(r,z) 


0 


L)(- 
2  MP, 


X)2 


+  (-T  +  T 
x  e 


J_wjl  °8n2  Tzr  ,-1/2 

_2Mp  ’  p  >  +  2.2J 

Z  *c  Kc  p  S 
rc 


(6.2) 
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In  the  filament  wound  composite,  the  radial  yield  stress  X  is 
essentially  that  of  the  resin  and  is  small  compared  with  the  axial  and 
circumferential  yield  stresses.  This  enables  the  approximation 


(i/e2)  <  <  (1/X2) 


(1/Z2)  <  <  (i/x2) 


(6.  3a,b) 


p  (r,z)  =  X (p  /p  )/[A)(— )  +  ~  (— )2  -  ^(^+  ^)],/2 
Koy  7  '  F°y  *V  'pc  Pc  2  SPC  pc  pc  pc 


(6-3c) 


From  either  of  Eqs.  (6.2)  or  (6.3)  it  can  be  clearly  seen  that  the 
transverse  shear  term  can  serve  to  decrease  the  yield  pressure. 
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7.  NUMERICAL  COMPUTATIONS  AND  DISCUSSION 

The  approximate  solution  to  the  three-dimensional  Navier  equa¬ 
tions  of  an  orthotropic,  infinite  circular  cylinder  requires  numerical 
verification  before  it  can  be  used.  For  this  reason  the  first  set  of 
numerical  computations  (which  were  all  performed  on  the  IBM  360/50  com¬ 
puter  located  at  the  Polytechnic  Institute  of  Brooklyn)  were  devoted  to 
comparing  the  approximate  results  developed  here  with  the  exact  solution 
to  a  transversely  isotropic  cylinder  subjected  only  to  periodically  spaced 
band  loads.  The  exact  results  were  given  to  the  author  by  H.  Levine,  who 
used  the  exact  analysis  described  in  Ref.  4.  The  cylinder  constants 
are 

Er-Efl-2E*  ' 

”  =  A  =  0.2 

(7-0 


Grz=Er/5.28  ,  Vzr=0.l6  ,  vr0=O.3O 


=  0.2 


1 


=  0.8 


6j  =  0.06l6 


The  comparison  is  shown  in  Figures  2  and  3.  The  results  indicate 
that  although  a  correction  term  [of  the  type  obtained  from  the  solution  of 
either  of  Eqs.  (4.15)  and  (4. 1 6) ]  is  desirable,  the  zeroth-approximation 
solution  gives  results  which  would  be  satisfactory  to  the  designer.  This 
solution  gives  a  fairly  accurate  representation  of  the  stress  distributions 
through  the  thickness  in  that  it  demonstrates  the  nonlinear  variation  of 
the  normal  stresses.  The  tendency  of  the  shear  stress  to  "peak"  in  the 
vicinity  of  the  load  discontinuity  is  predicted. 


31- 


The  second  set  of  computations  utilizes  the  approximate  three- 
dimensional  elasticity  solution  of  an  orthotropic  shell,  subjected  to  ex¬ 
ternal  pressure  and  an  axial  force  in  addition  to  prescribed  periodically 
spaced  band  loads.  The  results  offered  by  the  shell  theory  developed  in  Ref.  1 
are  compared  to  the  more  exact  results  obtained  here. 

In  Ref.  1,  the  band  loads  were  due  to  periodically  spaced  elastic 
ring  supports.  The  radial  displacement  of  a  ring  was  found  as  d  function 
of  the  magnitude  of  the  band  load  (pc)  in  a  separate  analysis  (either  ring 
theory  or  an  orthotropic  Lame  analysis).  The  ratio  PQ/PC  was  then  obtained 
by  matching  the  radial  displacements  of  the  ring  and  the  shell  at  the  ring- 
shell  interface.  In  the  present  analysis,  the  ratio  PD/PC  is  assumed  to 
have  the  same  numerical  values  as  those  determined  in  Ref.  1. 

II 

When  Flugge-type  shell  theory  was  applied  to  a  corresponding 

ring-supported  orthotropic  shell  it  was  found  that  p  /p  =  1.11.  The  addi- 

o  c 

tion  of  transverse  shear  deformation  to  this  analysis  resulted  in  p  /p  =  1.03* 

o  c 

Obviously,  the  results  of  the  present  three-dimensional  elasticity  solution  do  not 
reflect  the  restraining  effects  of  the  ring.  In  order  to  correct  this  it 
would  be  necessary  to  include  an  analysis  of  the  ring  and  then  match  dis¬ 
placements  at  the  ring-shell  interface. 

The  cylinder  constants  are 

6  £ 

E  =  2.49  x  1  O  ps i  ,  E  =6.14x10  psi  ,  E  =  4.74  x  10  psi 

r  9  z 

Grz  =  0.7  x  106psi  ,  vz0  =  0.136  ,  vQz  =  0.176 


Letter  dated  March  24,  1964  from  Mr.  W.P.  Couch,  David  Taylor  Model  Basin 


$£*  /  "?>  J'v***t^  *•*+*+& *#**  '■-  - 

■ftr 


£*  0.1755  ,  ~=  0.2417  ,  rQ  =  3-388  in.  ,  r,  =  3  in. 

6  j  =  0.0769 

X  =  0.22  x  I05psi  ,  0  =  1.7  x  lO^psi  ,  Z  =  1.0  x  JO^psi 


S  =  0.09  x  l(Tpsi 


(7-2) 


The  transverse  constants  £  ,  G  and  v  were  not  available,  as  they  are 
extremely  difficult  to  obtain  (Ref.  6)  and  were  simply  given  values  that 
were  felt  to  be  representative  of  this  type  of  material  (e.g.,  see  Ref.  7) • 

The  results  are  presented  in  Figures  4  to  10.  Shell  theory  gives 
excellent  results  at  midbay,  but  can  only  approximate  the  actual  state  of 
stress  at  the  frame.  The  yield  pressures  predicted  by  shell  theory  [using 
the  criterion  in  Eqs.  (6.1)  wi th  f  =  1]  are  fairly  accurate.  Although 
the  present  results  given  more  accurate  stress  distributions  for  the  band 
load  problem,  it  does  not  necessarily  follow  that  the  band  load  problem 
exactly  reflects  a  ring-shell  interaction  problem. 

Klosner  and  Levine  (Ref.  4)  found  that  transverse-shear-deformation 
shell  theory  did  not  lead  to  improvement  over  classical  shell  theory.  Such 
observations  can  also  be  made  from  the  present  calculations  (see  Figures 
4  and  5)-  Furthermore,  these  results  indicate  that,  at  a  sufficient  dis¬ 
tance  away  from  the  band  load,  the  shear  stress  can  be  approximated  with  a 
parabolic  function.  However,  the  results  displayed  in  Figs.  4  and  8  show 
that  the  axial  displacement  and  axial  stress  vary  cubical ly  (at  least) 
through  the  shell  thickness.  Hence,  in  the  vicinity  of  the  band  load,  where 
the  axial  displacement  reaches  its  largest  value,  plane  sections  do  not  re¬ 
main  plane  and  the  transverse  shear  stress  can  not  be  assumed  to  vary  para- 


>.¥  %  *«*!*.• 


33- 


bolically  through  the  thickness.  The  cubic  and  higher  order  thickness  terms 
combine  to  result  in  the  "peaking"  shown  in  Figure  6. 

The  Hill  yield  criterion  Eq.  (6.2)  predicts  that  yielding  begins 
at  a  fairly  low  pressure  (3500  psi)  near  the  inner  surface  at  the  load  dis¬ 
continuity.  Actual  tests  (Ref.  2)  indicate  that  the  shells  do  fail  at  the 
frame  near  the  inner  surface;  but  at  a  significantly  higher  pressure  (12200  psi). 
The  low  theoretical  yield  pressure  results  from  the  low  value  of  the  resin 
yield  stress.  When  this  yield  stress  was  increased  (i.e.,  X  was  set  equal 
to  1.0  x  105  psi)  the  lowest  value  of  the  theoretical  yield  pressure  was 
found  to  be  poy(rj,  L/2)  =  1 6323  ps i . 

One  might  speculate  that  the  resin  yields  at  a  very  low  pressure 
after  which  the  load  is  resisted  by  the  glass  fibers.  If  this  is  the  case, 
the  work  done  here  might  be  extended  to  a  layer  analysis  (see  Appendix), 
which  could  be  used  to  approximate  the  stresses  after  yielding  has  begun. 

It  appears  from  these  results  that  the  transverse  shear  stress 
distribution  has  very  little  effect  on  the  low  pressure  at  which  the  resin 
begins  to  yield  (which  was  also  predicted  by  the  shell  theory).  However, 
as  is  pointed  out  in  Ref.  1,  the  actual  stresses  in  the  individual  con¬ 
stituents  of  the  nonhomogeneous  filament  wound  composite  shells  must  be 
obtained  by  multiplying  the  stresses  obtained  here  by  suitable  stress- 
concentration  factors.  This  procedure  gives  rise  to  substantial  changes 
in  the  magnitude  of  the  initial  yield  pressure  (e.g.,  see  Ref.  1). 

It  is  interesting  to  note  that  the  shell  wall  gets  thicker 
under  pressure  (see  Fig.  9)-  This  can  be  traced  to  the  Poisson  effect 
of  the  rather  large  compressive  axial  and  circumferential  stresses. 


i , 

In  genera),  the  numerical  results  obtained  here  indicate  that  the 
approximate  elasticity  solution  developed  here  to  satisfy  the  orthotropic 
Navier  equations  gives  a  satisfactory  description  of  the  stress  distribution 
through  the  shell  thickness.  However,  before  this  technique  can  be  success¬ 
fully  applied,  the  transverse  elastic  constants  must  be  experimentally  deter¬ 
mined.  The  low  yield  pressures  obtained  here  may  correspond  to  the  "initial 
yielding"  state  described  by  Tsai  (Ref.  8)  in  which  the  load  deformation 
curve  of  a  similar  shell  subjected  to  internal  pressure  was  observed  to 
have  a  sudden  change  in  slope  at  a  low  pressure.  If  this  is  the  case,  the 
layer  analysis  suggested  in  the  Appendix  of  this  work  would  certainly  be 
appropriate. 


35. 


REFERENCES 


1.  Kempner,  Joseph;  Misovec,  A.P.  and  Herzner,  F.C.:  Ring-Stiffened 
Orthotropic  Circular  Cylindrical  Shell  under  Hydrostatic  Pressure. 

PIBAL  Rep.  No.  68-10,  Hay  1968. 

2.  Horn,  K.j  Buhl,  J.E.  and  Couch,  W.P.:  Hydrostatic  Pressure  Tests  of 
Unstiffened  and  Ring  Stiffened  Cylindrical  Shells  Fabricated  of  Glass 
Filament  Reinforced  Plastics,  David  Taylor  Model  Basin  Report  1745, 

Sept.  1963. 

3*  Klosner,  J.M.  and  Levine,  H.S.:  Further  Comparisons  of  Elasticity  and 
Shell  Theory  Solutions,  Polytechnic  Institute  of  Brooklyn,  PIBAL  Report 
No.  689,  July  1964. 

4.  Levine,  H.S.  and  Klosner,  J.M.:  Transversa 1 ly  Isotropic  Cylinders  under 
Band  Loads,  Jour.  Eng.  Mech.  Div.,  Proceedings  of  the  ASCE,  June  1967> 

P-  157- 

5-  Lekhnitskii,  S.G.:  Theory  of  Elasticity  of  an  Anisotropic  Elastic 
Body,  Holden  Day,  Inc.,  San  Francisco,  19&3- 

6.  Horn,  K.j  Couch,  W.P.  and  Wi liner,  A.R.:  Elastic  Macerial  Constants 
of  Filament  Wound  Cylinders  Fabricated  from  E-HTS/E787  and  S-HTS/E787 
Prepreg  Rovings,  David  Taylor  Model  Basin  Report  1823,  Feb.  1 966. 

7.  Myers,  N.C.j  Lee,  G.D.j  Wright,  F.C.  and  Daines,  J.V.:  Investigation 
of  Structural  Problems  with  Filament  Wound  Deep  Submersibles,  Final 
Report,  H.l.  Thompson  Fiber  Glass  Co.,  January  1964. 

8.  Tsai,  S.W.j  Adams,  D.F.  and  Doner,  D.R.:  Analyses  of  Composite 
Structures,  NASA  CR  620,  November  1966. 

9.  Hill,  R.j  A  Theory  of  the  Yielding  and  Plastic  Flow  of  Anisotropic 
Metals,  Proceedings  Royal  Society,  A 1 93,  pp.  281-286,  1948. 


9* 


36. 


•  •frxagvz?:  ■ 


APPENDIX 


!t  may  be  desirable  to  prescribe  a  variety  of  combinations  of 
radial  stress,  shear  stress,  axial  displacement  and  radial  displacement  as 
periodic  functions  of  z  on  either  of  the  shell  surfaces  (for  example,  the 
ring  supported  shell  may  be  analysed  by  matching  ring  and  shell  displacements) 
Therefore  a  more  general  discussion  of  possible  boundary  conditions  is  offered 
On  each  shell  surface  (outer  and  inner),  any  two  of  the  following 


four  quantities  must  be  prescribed: 

4  _xo 

e  [«  '(VVi2>  ’  ai}HiilG1iEn  =  H" 

1  =  1  2  4 


(radial  stress) 


(Ala) 


e  («  Vr.-'lCn'ii  »h"2 

1  =  1  2  J 


n 

U"  =  G 


E  H..G"  =  Fn 

M  "  11 


(transverse  shear  stress) 


(radial  displacement) 


(Alb) 


(Ale) 


(axial  displacement)  (Aid) 


where  H?,  Gn  and  Fn  are  Fourier  series  coefficients  obtained  by  expanding 
the  prescribed  functions  of  z. 


It  is  possible  to  perform  a  layer  analysis  by  matching  the  four 
quantities  listed  in  Eqs.  (Al)  across  each  interface.  If  a  shell  consists 
of  m  layers,  4m  arbitrary  constants  would  have  to  be  determined  The  inter¬ 
action  equations,  which  would  be  4m~4  in  number,  could  be  obtained  from 
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(H")  = 

j-» 

(Hj). 

(A2a) 

i 

J 

(H2> .  .  " 

(HP*) 

✓ 

(A2b) 

J-t 

J 

Gn  ,  =  Gn 
J-i  J 

(A2c) 

F°  ,  =  Fn 
J-'  J 

(A2d) 

where  j  =  2,  **',  m.  The  subscripts  1  and  m  correspond  to  the  inner  and 
outer  shell  layers,  respectively.  The  remaining  4  equations  are  found  by 
applying  two  of  Eqs .  (A!)  on  layers  1  and  m.  The  corresponding  generalized 
plane  strain  problem  must  also  be  solved  layer  by  layer  with  the  radial 
stress  and  both  displacements  being  matched  at  each  interface. 

If  the  layers  were  permitted  an  axial  motion  relative  to  each 
other  (as  might  occur  after  initial  yielding)  Eqs.  (A2b)  and  (<\2d)  could  be 
replaced  by 


(H")  =  Tn 

(A2'b) 

j-1 

(H")  =  Tn 

(A2'd) 

where  Tn  would  be  a  constant  corresponding  to  a  maximum  shear  stress  at 


resin  yield. 
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SHELL  GEOMETRY  AND  LOADING 
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FIG.  2  SOLUTIONS  OF  A  TRANSVERSALLY  ISOTROPIC  CYLINDER 

(AXIAL  AND  TRANSVERSE  SHEAR  stress) 
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FIG.  7 

RADIAL  STRESSES  IN  ORTHOTROPIC  CYLINDER 
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FIG.  8  AXIAL  DISPLACEMENTS  OF 
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FIG.  9  RADIAL  DISPLACEMENTS  OF 
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